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Abstract A number of different approaches to quantum gravity are at least partly 
phenomenologically characterized by their treatment of Lorentz symmetry, in par- 
ticular whether the symmetry is exact or modified/broken at the smallest scales. For 
example, string theory generally preserves Lorentz symmetry while analog gravity 
and Lifshitz models break it at microscopic scales. In models with broken Lorentz 
symmetry there are a vast number of constraints on departures from Lorentz invari- 
ance that can be established with low energy experiments by employing the tech- 
niques of effective field theory in both the matter and gravitational sectors. We shall 
review here the low energy effective field theory approach to Lorentz breaking in 
these sectors and present various constraints provided by available observations. 



1 Introduction 

Our understanding of the observed laws of Nature is currently based on two different 
theories: the Standard Model of Particle Physics (SM), and General Relativity (GR). 
However, in spite of their phenomenological successes, SM and GR leave many 
fundamental theoretical questions unanswered. First of all, part of the success of 
the SM has been the recognition that symmetry breaking is an important part of 
modern physics and that what appear to be multiple forces at low energies can often 
be described in a unified manner. The prime example of this is the Glashow-Salam- 
Weinberg theory of electroweak interactions. Since such unification is possible, and 
since many physicists feel that our understanding of the fundamental laws of Nature 
is deeper and more accomplished if we are able to reduce the number of degrees of 
freedom in a theory, much effort has been spent trying to construct unified theories 
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in which not only all the sub-nuclear forces are seen as different aspects of a unique 
interaction, but also gravity is included in a consistent manner as merely part of the 
overall structure. 

Another important reason why we seek for a new theory of gravity comes di- 
rectly from the gravity side. We know that GR fails to be a predictive theory in 
some regimes. Indeed, many solutions of Einstein's equations are singular in some 
region, and GR is not able to make any prediction in those regions of spacetime. 
Moreover, there are honest classical solutions of the Einstein equations that contain 
closed time-like curves, which would allow traveling back and forth in time with the 
associated causal paradoxes. Finally, the problem of black-hole evaporation con- 
sidered just within the framework of semi-classical gravity clashes with quantum 
mechanical unitary evolution. 

This long list of puzzles spurred intense research toward a quantum theory of 
gravity that started almost immediately after Einstein's proposal of GR and which 
is still one of the most active areas of theoretical physics. The quantum gravity 
problem is not only conceptually and technially challenging, it has also been an 
almost metaphysical pursuit for several decades, in that what progress has been 
made is on the theoretical side and the experimental aspect has been (for good 
reason) neglected. Indeed, we expect QG effects at experimentally/observationally 
accessible energies to be extremely small, due to suppression by the Planck scale 
M p i = ^Tic/Gn — 1-22 x 10 19 GeV/c 2 . In this sense it has been considered (and it 
is still considered by many) that only ultra-high-precision (or Planck scale energy) 
experiments would be able to test quantum gravity models. 

It was however realized (mainly over the course of the past decade) that the 
situation is not quite as bleak as it appears. In fact, quantum gravitational models 
beyond GR have shown that there can be several of what we term low energy "relic 
signatures" of quantum gravitational effects which would lead to deviation from the 
standard theory predictions (SM plus GR) in specific regimes. Some of these new 
phenomena, which comprise what is often termed "QG phenomenology", include: 

• Quantum decoherence and state collapse [ 1 ] 

• QG imprint on initial cosmological perturbations [2] 

• Cosmological variation of couplings [3, 4] 

• TeV Black Holes, related to extra-dimensions [5] 

• Violation of discrete symmetries [6] 

• Violation of space-time symmetries [7] 

In this chapter we will focus upon the phenomenology of violations of space-time 
symmetries, and in particular of Local Lorentz invariance (LLI), a pillar both of 
quantum field theory as well as GR (LLI is a crucial part of the Einstein Equivalence 
Principle on which metric theories of gravity are based). 
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Contrary to the common perception, explorations of the possible breakdown of LLI 
have a long standing history. It is however undeniable that the last twenty years 
have witnessed a striking acceleration in the development both of theoretical ideas 
as well as of phenomenological tests previously unimagined. We shall here present 
an admittedly incomplete review of these developments. 



2.1 Early works 

The possibility that Lorentz invariance violation (LV) could play a role in physics 
dates back at least sixty years [8, 9, 10, 11, 12, 13] and in the seventies and eighties 
there was already a well established literature investigating the possible phenomeno- 
logical consequences of LV (see e.g. [14, 15, 16, 17, 18, 19]). 

The relative scarcity of these studies in the field was due to the general expec- 
tation that new effects would only appear in particle interactions where the particle 
energies were of order the Planck scale, as only at those energies would the natural 
QG suppression by powers of the Planck scale be overcome. However, it was only 
in the nineties that it was clearly realized that there are special situations in which 
new effects, even if highly suppressed, can have observational consequences. These 
situations were termed "Windows on Quantum Gravity". 



2.2 The dawn of Quantum Gravity Phenomenology 

At first glance, it appears hopeless to search for effects suppressed by the Planck 
scale. Even the most energetic particles ever detected (Ultra High Energy Cosmic 
Rays, see, e.g., [20, 21]) have E < 10 11 GeV ~ 10~ 8 M p /. However, tiny corrections 
can be magnified into an observational effect if the physics in question involves 
not just the Planck scale and the energy scale of the particle, but also another scale 
such as the mass of a light particle or a cosmological travel time. In these situations 
observables can be constructed that leverage the scales against one another to create 
an actually measurable effect at achievable energies (i.e. anything from the energies 
in tabletop optical experiments to cosmic rays. See, e.g., [7, 22] for an extensive 
review). 

A partial list of these windows on QG includes: 

• sidereal signal variations as a laboratory apparatus such as an optical cavity 
moves with respect to a preferred frame or direction 

• cumulative effects: long baseline dispersion and vacuum birefringence (e.g. of 
signals from gamma ray bursts, active galactic nuclei, pulsars) 
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• anomalous (normally forbidden) threshold reactions allowed by LV terms (e.g. pho- 
ton decay, vacuum Cerenkov effect) 

• shifting of existing threshold reactions (e.g. photon annihilation from Blazars, 
ultra high energy protons pion production) 

• LV induced decays not characterized by a threshold (e.g. decay of a particle from 
one helicity to the other or photon splitting) 

• maximum attainable particle velocities different from c (e.g. synchrotron peak 
from supernova remnants) 

• dynamical effects of LV background fields (e.g. gravitational coupling and addi- 
tional wave modes) 

It is rare one can assign a definitive "paternity" to a field, and our so called 
"Quantum Gravity Phenomenology" is no exception. However, among the papers 
commonly accepted as seminal we can cite the one by Kostelecky and Samuel [23] 
that already in 1989 envisaged, within a string field theory framework, the pos- 
sibility of non-zero vacuum expectation values (VEV) for some Lorentz breaking 
operators. This work led later on to the development of a systematic extension of 
the SM (what was later on called "minimal standard model extension" (mSME)) 
incorporating all possible Lorentz breaking, power counting renormalizable opera- 
tors (i.e. of mass dimension < 4) by Colladay and Kostelecky [24]. This provided 
a framework for computing in effective field theory the observable consequences 
for many experiments and led to much experimental work setting limits on the LV 
parameters in the Lagrangian (see e.g. [25] for a review). 

Another seminal paper was that of Amelino-Camelia and collaborators [26] 
which highlighted the possibility to cast observational constraints on Planck-suppressed 
violations of Lorentz invariance in the photon dispersion relation by examining the 
propagation of light from remote astrophysical sources like gamma ray bursters 
(GRBs) and active galactic nuclei (AGN). Finally, we also mention the influential 
papers by Coleman and Glashow [27, 28, 29] which brought the subject of sys- 
tematic tests of Lorentz violation to the attention of of the broader community of 
particle physicists. 

Let us stress that this is necessarily an incomplete account of the literature that 
investigated departures from Special Relativity. Several papers appeared in the same 
period and some of them anticipated many important results, see e.g. [30, 31]; unfor- 
tunately at the time of their appearance they were hardly noticed (and seen by many 
as too "exotic"). 

In the first decade after 2000 the field reached a concrete maturity and many 
papers pursued both a systematization of the various frameworks and the available 
constraints (see e.g. [32, 33, 34]). In this sense another crucial contribution was the 
development of an effective field theory approach also for higher order (mass di- 
mension greater than four), naively non-power counting renormalizable, operators 
l . This was first done for rotationally invariant dimension 5 operators in QED [38] 
by Myers and Pospelov which was later on extended to larger sections of the stan- 



1 Anisotropic scaling [35, 36, 37] techniques were recently recognized to be the most appropriate 
way of handling higher order operators in Lorentz breaking theories and in this case the highest 
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dard model by Bolokhov and Pospelov [39] and dimension 6 operators by Mattingly 
[40] . The general set of higher dimension operators for free photons and fermions 
have recently been cataloged by Kostelecky and Mewes [41, 42]. 

Why did all this attention to Lorentz breaking frameworks and observations de- 
velop in the late nineties and in the first decade of the new century? The answer 
is twofold as it is related to important developments coming from experiments and 
observation as well as from theoretical investigations. Observationally, there were 
a number of puzzling observations related to gravity that spawned a correspond- 
ing growth in the zoo of quantum gravity models/scenarios with a low energy phe- 
nomenology. For example, in cosmology these are the years of the striking realiza- 
tion that our universe is undergoing an accelerated expansion phase [43, 44] which 
apparently requires a new exotic cosmological fluid, called dark energy, which vi- 
olates the strong energy condition (to be added to the already well known, and still 
mysterious, dark matter component). 

Also in the same period high energy astrophysics provided some new observa- 
tional puzzles directly related to Lorentz symmetry, first with the apparent absence 
of the Greisen-Zatsepin-Kuzmin (GZK) cut off [45, 46] (a suppression of the high- 
energy tail of the UHECR spectrum due to UHECR interaction with CMB photons) 
as claimed by the Japanese experiment AGASA [47], and later on with the so called 
TeV-gamma rays crisis, i.e. the apparent detection of a reduced absorption of TeV 
gamma rays emitted by AGN [48]. Both these "crises" later on subsided or at least 
alternative, more orthodox explanations for them were advanced. However, their 
existence undoubtedly boosted the research in the field at that time. 

It is perhaps this past "training" that made several collaborations within the quan- 
tum gravity phenomenology community strongly emphasize the apparent incom- 
patibility of the recent CERN-LNGS based experiment OPERA [49] measurement 
of superluminal propagation of muonic neutrinos with Lorentz violating EFT (see 
e.g. [50, 51, 52, 53]. There is now evidence that the Opera measurement might be 
flawed due to unaccounted experimental errors and furthermore has been refuted by 
a similar measurement of the ICARUS collaboration [54]. Nonetheless, this claim 
propelled a new burst of activity in Lorentz breaking phenomenology which might 
still provide useful insights for future searches. 

Parallel to these exciting developments on the experimental/observational side, 
theoretical investigations provided new motivations for Lorentz breaking searches 
and constraints. Indeed, specific hints of LV arose from various approaches to Quan- 
tum Gravity. Among the many examples are the above mentioned string theory 
tensor VEVs [23] and space-time foam models [55, 26, 56, 57, 58], then semi- 
classical spin-network calculations in Loop QG [59], non-commutative geome- 
try [60, 61, 62], some brane-world backgrounds [63]. 

More recently, it cannot be omitted the role associated with the development of 
Lorentz breaking theories of gravity from early studies [64, 65, 66, 67, 68] to more 
systematic approaches such as Einstein- Aether [69, 70, 71] and Hofava-Lifshitz 
[36] gravity. Finally, there was the vigorous development over the same time of the 



order operators are indeed crucial in making the theory power counting renormalizable. This is 
why we shall adopt sometimes the expression "naively non renormalizable" 
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so called condensed matter analogues of "emergent gravity" [72], which showed 
how approximate local Lorentz invariance can arise from a fundamentally Galilean 
theory. 

Many of these approaches yield a low energy description of Lorentz violation 
in terms of effective field theory, and so before we delve into the specific opera- 
tors that people have considered we remark on some generic aspects of embedding 
phenomenologically acceptable Lorentz violation in a low energy EFT. 



3 Modified dispersion relations and their naturalness 

We will concentrate here on the free field modifications to EFT from LV and hence 
talk primarily about modifications to dispersion relations. This restriction isn't as 
limiting as it might seem. For example, the necessary LV interaction terms generated 
when one wants to maintain gauge covariance in a LV theory [24] are controlled 
by the same LV coefficients that control the free theory so no new LV coefficients 
are introduced. Since most experimental work is sensitive to free field behavior, 
constraints on those coefficients are generated from the phenomenology of the free 
part of the Lagrangian rather than the interaction part. One can of course add other 
LV interaction terms by hand - such terms have not been as extensively studied in 
the literature and so we will not focus on them here. There is one caveat: when one 
sees constraints in the literature from modified particle decay mechanisms there is 
usually an assumption that the interactions of the standard model hold with only 
small modifications if any. This is reasonable, as if there is only a tiny modification 
to the free field equations then by the argument above about how gauge generated 
LV interaction terms are controlled by the same coefficients as the free field part, 
the corresponding modifications to the interaction terms would also be small. This 
would, in general, modify the rate of various reactions but the rate difference given 
the constraints detailed below would be unobservable. 

Turning back to the issues of modified dispersion relations in LV EFT, one has 
to firstly recognize that calculations of a specific dispersion modification from a 
specific quantum gravity theory are in general problematic and that cases where one 
can do so [55] are the exception rather than the norm. A more reasonable approach is 
to therefore simply consider a generic momentum expansion of a dispersion relation 
in a specified observer's frame of 

N=l 

where the low energy speed of light c = 1 and each rj W is an arbitrary rank N tensor 
with mass dimension 2—N and the a, indices run over spacetime coordinates. 

Let us note that this ansatz assumes that the propagating mass eigenstates are 
also eigenstates of the Lorentz violating physics. This does not need to be the case, 
and having the eigenstates of Lorentz violation not match the mass eigenstates can 
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be useful when trying to analyze the effects of Lorentz violation on neutrinos. How- 
ever, since such a mismatch in other sectors would introduce oscillations and other 
unseen effects for particles other than neutrinos, i.e. those particles we are primarily 
concerned about. Therefore we shall not consider such possibilities outside of the 
neutrino sector. 

The jj^'s can be mapped on to coefficients in a corresponding Lagrangian, al- 
though we note for the reader that in a generic Lagrangian there are many other 
coefficients that don't influence free field dispersion (c.f. [41]). Needless to say, as 
N increases, there are a multitude of possible dispersion relations. Testing all such 
possible combinations can be done, of course, and would be the most systematic 
way to evaluate the possibility of Lorentz violation. 

One can, however, simplify the possible set of dispersion by making various as- 
sumptions such as CPT invariance, rotational invariance, etc. Rotational invariance 
is one of the most common assumptions made, for three reasons. 2 The first is that 
it dramatically reduces the number of possible dispersion terms while still preserv- 
ing interesting phenomenology. Second, in many instances rotational invariance is 
more stringently tested than boost invariance so it provides a good model for test- 
ing boost invariance alone. Finally, many QG models out there that do not predict 
exact Lorentz invariance still preserve rotational invariance. (See [34] for further 
discussion about this assumption.) 

With rotational invariance in mind in some frame, a common assumption for the 
dispersion relation is then 

E 2 =p 2 +m 2 + £ f[ N p N , (2) 

N=\ 

where p is the magnitude of the three momentum. This type of expansion assumes 
high energy particles and that the corrections are small, so anywhere E would have 
appeared on the RHS of (2) it is replaced by p. In general one can allow the LV 
parameters rjw to depend on the particle type, and indeed it turns out that they must 
sometimes be different but related in certain ways for photon polarization states, 
and for particle and antiparticle states, if the framework of effective field theory is 
adopted. The lowest order LV terms (p, p 2 , p 3 , p 4 ) have been the terms that have 
generated the most attention (c.f. [7, 73] and references therein) 3 . 



2 A notable exception to this assumption is the Standard Model Extension and associated tests. 
Rotational invariance is not assumed in the program of the standard model extension as it considers 
all terms at each mass dimension. 

3 We disregard here the possible appearance of dissipative terms [74] in the dispersion relation, 
as this would correspond to a theory with unitarity loss and to a more radical departure from 
standard physics than that envisaged in the framework discussed herein (albeit a priori such dissi- 
pative scenarios are logically consistent and even plausible within some quantum/emergent gravity 
frameworks). 
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3.1 The naturalness problem 

From a EFT point of view the only relevant operators should be the lowest order 
ones, i.e. those of mass dimension 3,4 corresponding to terms of order p and p 2 
in the dispersion relation. Situations in which higher order operators contribute to 
the dispersion as much as the lowest order ones at some energy are only possible 
at the cost of a severe, indeed arbitrary, fine tuning of the coefficients tjaj (which 
we discuss below). However, as we shall see current observational constraints are 
incredibly tight on dimension 3 operators and very severe on dimension 4 ones. This 
is kind of obvious, given that these operators would end up modifying the dispersion 
relation of elementary particles at easily achievable energies. Dimension 3 operators 
would dominate as p — > while the dimension 4 ones would generically induce a 
species dependent, constant shift in the limit speed for elementary particles. Hence 
one is left with two less than perfect approaches. 

First, one can assume the standard EFT hierarchy, stop testing at operators of 
mass dimension three and four, and, due to the tightness of limits, argue that Lorentz 
invariance is likely an exact symmetry of nature and that QG/emergent models that 
do not respect the symmetry should be discounted. Or, one can avoid assumptions 
about whether any additional new physics comes into play between everyday ener- 
gies and the QG scale and so not assume a particular hierarchy between operators of 
various mass dimension. This is the approach many phenomenologists take: simply 
start at the lowest dimension operators, derive constraints, and work upwards in N 
as far as one can observationally go without imposing any necessary hierarchy. 

Not assuming a necessary hierarchy (3) and simply constraining the coefficients 
T]n at each order is perfectly good phenomenologically, and we will take that ap- 
proach going forward, but it is important that the reader understand why it is the- 
oretically unnatural. The reason such a hierarchy is unnatural is simple: in EFT 
radiative corrections will generically allow the percolation of higher dimension 
Lorentz violating terms into the lower dimension terms due to the interactions of 
particles [75, 76]. 

In EFT loop integrals will be naturally cut-off at the EFT breaking scale, if such 
scale is as well the Lorentz breaking scale the two will effectively cancel leading 
to unsuppressed, coupling dependent, contributions to the base dimension four ki- 
netic terms that generate the usual propagators. Hence radiative corrections will not 
allow a dispersion relation with only p 3 or p A Lorentz breaking terms but will auto- 
matically induce extra unsuppressed LV terms in p and p 2 which will be naturally 
dominant. One could argue that RG effects might naturally suppress the sizes of 
these coefficients at low energies. As we shall see, in specific models where the RG 
flow has been calculated, the running of LV coefficients is only logarithmic and so 
there is no indication that RG flow will actually drive coefficients to zero quickly in 
the infrared. 

Several ideas have been advanced in order to resolve such a "naturalness prob- 
lem" (see eg. [34]). While it would be cumbersome to review all the proposals here, 
we point out two of the more prominent ideas, both of which involve introducing 
a new scale into the problem in addition to the Lorentz violating scale. If there are 
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two scales M and ;ii involved, there can be a hierarchy of LV coefficients different 
than the naive one, for example 

U 2 u 1 

(i n M , u u M , n u M V J 

or ^ 

m = m^, ^4 = 774^ (4) 

where M is the QG scale, usually taken to be the Planck scale, jj, is some other far 
lower energy scale, and tja? is a de-demensionalized coefficient usually assumed to 
be O(l). The exact hierarchy, whether it involves terms of every mass dimension as 
in (3), or only even dimension as in (4) (which can be accomplished up to = 4 
by imposing CPT) is model dependent. We now briefly describe two ideas that have 
been put forward that would generate such scales. 



3.1.1 A new symmetry 

Most of the aforementioned proposals implicitly assume that the Lorentz breaking 
scale is the Planck scale. One then needs the EFT scale (which can be naively iden- 
tified with what we called previously jj.) to be different from the Planck scale and 
actually sufficiently small so that the lowest order "induced" coefficients can be sup- 
pressed by suitable small rations of the kind fi p /M q where p,q are some positive 
powers. 

A possible solution in this direction can be provided by introducing what is com- 
monly called a "custodial symmetry" - a symmetry other than Lorentz that forbids 
lower dimension operators and is broken at the low scale jj,. The most plausible 
candidate for this role is supersymmetry (SUSY) [77, 78]. SUSY is by definition 
a symmetry relating fermions to bosons i.e. matter with interaction carriers. As a 
matter of fact, SUSY is intimately related to Lorentz invariance. Indeed, it can be 
shown that the composition of at least two SUSY transformations induces space- 
time translations. However, SUSY can still be an exact symmetry even in presence 
of LV and can actually serve as a custodial symmetry preventing certain operators 
to appear in LV field theories. 

The effect of SUSY on LV is to prevent dimension < 4, renormalizable LV oper- 
ators to be present in the Lagrangian. Moreover, it has been demonstrated [77, 78] 
that the renormalization group equations for Supersymmetric QED plus the addi- 
tion of dimension 5 LV operators a la Myers & Pospelov [38] do not generate lower 
dimensional operators if SUSY is unbroken. However, this is not the case for our 
low energy world, of which SUSY is definitely not a symmetry. 

The effect of soft SUSY breaking was also investigated in [77, 78]. It was found 
there that, as expected, when SUSY is broken the renormalizable operators are gen- 
erated. In particular, dimension K ones arise from the percolation of dimension K + 2 
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LV operators . The effect of SUSY soft-breaking is, however, to introduce a sup- 
pression of order m 2 jMp\ (k = 3) or (m s /Mpi) 2 (K = 4), where m s ~ 1 TeV is the 
scale of SUSY soft breaking. Although, given present constraints, the theory with 
K = 3 needs a lot of fine tuning to be viable, since the SUSY-breaking-induced sup- 
pression is not enough powerful to kill linear modifications in the dispersion relation 
of electrons, if K = 4 then the induced dimension 4 terms are suppressed enough, 
provided m s < 100 TeV. Current lower bounds from the Large Hadron Collider are 
at most around 950 GeV for the most simple models of SUSY [79] (the so called 
"constrained minimal supersymmetric standard model", CMSSM). 

Finally, it is also interesting to note that the analogue model of gravity can be used 
as a particular implementation of the above mentioned mechanism for avoiding the 
so called naturalness problem via a custodial symmetry. This was indeed the case of 
multi-BEC [80,81]. 



3.1.2 Gravitational confinement of Lorentz violation 

The alternative to an extra symmetry is to turn the problem upside down and posit 
that the Lorentz breaking scale (the M appearing in the above dispersion rela- 
tions) is not set by the Planck scale, but is instead the lower scale jj,. If one does 
this and begins with a theory which has higher order Lorentz violating operators 
only in the gravitational sector, then one can hope that the gravitational coupling 
Gn ~ Mpj 2 will suppress the "percolation" to the matter sector where the constraints 
are strongest. Matter Lorentz violating terms will all possess factors of the order 
(fx /Mpi) 2 which can become strong suppression factors if ji <C Mp\. This is the idea 
underlying the work presented in [82] which applies it to the special case of Hofava- 
Lifshitz gravity. There it was shown that indeed a workable low energy limit of the 
theory can be derived through this mechanism which apparently is fully compati- 
ble with existing constraints on Lorentz breaking operators in the matter sector. In 
our opinion, this new route deserves further attention and should be more deeply 
explored in the future. 



4 Dynamical Frameworks I: Rotationally Invariant CPT even 
Standard Model Extension with mass dimension five and six 
operators 

We now turn from theoretical considerations about naturalness and retreat to the 
phenomenological approach of "constrain everything as best one can observation- 
ally". There are various systematic frameworks and approaches for this. The "Stan- 
dard Model Extension" or SME contains all possible Lorentz violating tensors that 
can be coupled to standard model fields without changing the field content or vio- 



4 We consider here only K" — 3,4, for which these relationships have been demonstrated. 
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lating gauge symmetry. The standard model extension can be split into the "Min- 
imal Standard Model Extension" [23] of Kostelecky et. al., which contains only 
renormalizable operators, and the full SME, which contains the infinite tower of 
non-renormalizable higher mass dimension operators. As one can imagine, there 
are dozens of possible operators even for the renormalizable case, and the number 
of operators up to mass dimension six is in the hundreds. Many of the operators can 
be constrained by similar methods, so for the purposes of this introduction we will 
concentrate on the most well studied sector of the SME, that of rotationally invari- 
ant QED. In particular we will concentrate on the interactions of photons, electrons, 
and protons. 

Usually when defining a field theory one starts with the renormalizable operators 
and proceeds in increasing mass dimension. Here we start with higher mass dimen- 
sion and work backwards for purely pedagogical reasons - many of the reactions 
and constraints we describe in detail for the higher dimension operators will also be 
useful for setting constraints on various lower dimension operators. The first mass 
dimension we explore is mass dimension six, as this is the highest mass dimension 
where all the operators have been classified and significantly studied. We break the 
operators into CPT even and odd classes as a different set of observations can be 
used to constraint CPT odd operators. 



4.1 The model 

A list of CPT even, rotationally invariant mass dimension five and six LV terms was 
computed in [40] through the same procedure used by Myers & Pospelov for di- 
mension 5 LV (see below), and this has been extended to non-rotationally invariant 
operators in [41, 42]. With rotation invariance all LV tensors must reduce to suitable 
products of a time-like vector field, usually denoted u a . This is usually taken to be 
unit, so that in the frame of the observer whose world line is tangent to u a , u a has 
components (1,0,0,0). This allows us to express constraints solely in terms of the 
numerical coefficient involved in any M a -matter interaction term. Of course, the ac- 
tual direction of u a is technically arbitrary. However, the common choice, which we 
make here, is to define u a to be aligned with the rest frame of the cosmic microwave 
background. In terms of u a the known mass dimension six fermion operators are 

--L.y(u-D) 3 (u ■ y)(a| 6) P L + a^P R ) ¥ (5) 

_ ' ¥(m . D )u(u ■ y){a^P L + af ] P R )y , 
M n 

where Pr± are the usual left and right spin projectors Pr.l = (1 ± 7 5 )/2 and D is 
the gauge covariant derivative. All the coefficients a are dimensionless because we 
factorize suitable powers of the Planck mass out explicitly. In addition there is a 
CPT even dimension five term [40], 
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^y(u-D) 2 {a { L 5) P L + a { R 5) P R )\i/. (6) 



The known photon operator is 

-^ , F^u^u a {u-d)F av . (7) 



1 nWtfv,, 
2M 2 



j pi 



From these operators, the dispersion relations of fermions and photons can be 
computed, yielding 



(6)r-3 (6) -3 

E 2 —p 2 — m 2 = ^(e + s P ) + (e - s P ) + £- (4 5) + 4 y + 4 5) «| 5) 

M pl A/ pl Mpi 



G) 2_ yt 2 = j3 (6) 

M, 2 



A- 4 
PI 



where m is the electron mass and where s = a ■ p/|p| is the helicity of the fermions. 
The a terms contribute as m 2 /Mp,, i.e. highly suppressed, and so will be neglected. 

In general, a LV dispersion for a particle with a certain set of quantum numbers 
(mass, spin, etc.) will be of the form E 2 = p 2 + /W p" /M^ 2 , and so we will often 
refer to type "n" LV. For example, because the high-energy fermion states are almost 
exactly chiral, we can further simplify the fermion dispersion relation eq. (8) (with 
R = +,! = -) 

E 2 =p 2 + m 2 + — T,(V + !ji 4) -4- W 
Mp\ (± Mp, 

where i]W is the dispersion coefficient of the LV p" term in the dispersion relation 
for the fermion. We choose rj as the coefficient as this nomenclature is common 
in the literature. Similarly, £W will refer to the generic dispersion coefficient for a 
photon (so in the case above ^' 4 ' = /3' 6 '). As it is suppressed by a factor of order 
m/Mpi, we will drop the quadratic modification generated by the dimension five 
operator. Indeed this can be safely neglected, provided that E > \/mMp\. Let us 
stress however, that this is exactly an example of a dimension 4 LV term with a 
natural suppression, which for electron is of order m e /Mp\ ~ 10~ 22 . Therefore, any 
limit larger than 10~ 22 placed on this term would not have to be considered as an 
effective constraint (to date, the best constraint for a rotationally invariant electron 
LV term of dimension 4 is <9(1CT 16 ) [83]). Note that modulo this the CPT even 
dimension five operator for fermions has the same effect on the dispersion as the 
CPT even dimension six in that it generates a p 4 term, so we will generally just write 
constraints directly on T]± . It may seem puzzling that in a CPT invariant theory we 
distinguish between different fermion helicities in (9). However, although they are 
CPT invariant, some of the LV terms displayed in eq. (7) are odd under P and T. 

CPT invariance allows us to determine a relationship between the LV coefficients 
of fermions and anti-fermions. Indeed, to obtain these we simply realize that, by 
CPT, the dispersion relation of the anti-fermion is given by (8), with the replace- 
ments s —> —s and p —p. If q,q denote a charge fermion and anti-fermion, then 
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the relevant anti-fermion coefficient rj^ is such that T|g = 7]^. , where q ± indi- 
cates an anti-fermion of positive/negative helicity (and similarly for the q±). Let us 
anticipate that the same argument used above leads to the conclusions that for disper- 
sion relations with odd powers of "n" (e.g. p 3 type dispersion relations) one obtains 
^Krfrf) _ _j^('Wrf^ jj ence f or arbitrary "n" one would expect rji*' = (-1)"tj^"'. 
This different behavior between even and odd powers "n" type dispersion relations 
leads to quite distinct phenomenologies as we shall see later. 



4.2 Constraints 

4.2.1 Threshold reactions 

Threshold reactions with ultra-high energy cosmic rays provide the only significant 
constraints on the above operators. A threshold reaction is a reaction that does not 
occur above a certain energy scale, which we call the "threshold energy". An in- 
teresting and useful phenomenology of threshold reactions is introduced by LV in 
EFT; also, threshold theorems can be rederived [33]. Sticking to the present case 
of rotational invariance and monotonic dispersion relations (see [84] for a general- 
ization to more complex situations), the kinematics of threshold reactions yield a 
number of useful phenomenological facts about these reactions [32]: 

• Threshold configurations still correspond to head-on incoming particles and par- 
allel outgoing ones 

• The threshold energy of existing threshold reactions can shift, and upper thresh- 
olds (i.e. maximal incoming momenta at which the reaction can happen in any 
configuration) can appear 

• Pair production can occur with unequal outgoing momenta 

• New, normally forbidden reactions can be viable 

LV corrections are surprisingly important in threshold reactions because the LV 
term (which as a first approximation can be considered as an additional energy de- 
pendent "mass") should be compared not to the momentum of the involved particles, 
but rather to the (invariant) mass of the heaviest particle in the interaction. Thus, an 
estimate for the threshold energy is 

( m 2 M n-2\ >/« 

^pi^r) • (10) 

where m is the mass of the heaviest particle involved in the reaction. Interesting 
values for p t h are discussed, e.g., in [32] and given in Tab. 1. Reactions involving 
neutrinos are the best candidate for observation of LV effects, whereas electrons and 
positrons can provide results for n = 3 theories but cannot readily be accelerated by 
astrophysical objects up to the required energy for n = 4. In this case reactions of 
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Table 1 Values of p^, according to eq. (10), for different particles involved in the reaction: neu- 
trinos, electrons and proton. Here we assume 7]'"' ~ 1. 





m v -0.1 eV 


m e ~ 0.5 MeV 


m p ~ 1 GeV 


n = 2 


0.1 eV 


0.5 MeV 


1 GeV 


n = 3 


500 MeV 


14 TeV 


2PeV 


n = 4 


33 TeV 


74PeV 


3 EeV 



protons can be very effective, because cosmic rays can have measured energies well 
above 3 EeV. We now discuss two threshold reactions are of particular use when 
constraining n = 4 LV. 

LV-allowed threshold reactions: y-decay 

The decay of a photon into an electron/positron pair is made possible by LV because 
energy-momentum conservation may now allow otherwise forbidden reactions to 
occur. Since the decay is a reaction described by the fundamental QED vertex, the 
rate once above threshold will be quite fast. The threshold for this process, is set by 
the condition [34] 

/ ro 2 Afgf 2 V 7 " 

where is a linear combination of tj Notably, the electron- 

positron pair can be created with slightly different outgoing momenta (asymmetric 
pair production). Furthermore, the decay rate is extremely fast above threshold [34] 
and is of the order of (10 ns) _1 (« = 3) or (10~ 6 ns) _1 (n = 4) if the LV coefficients 
are of O(l). 

LV-allowed threshold reactions: Vacuum Cerenkov 

In the presence of LV, the process of Vacuum Cerenkov (VC) radiation — > q^y, 
where q is a charged fermion can occur as this is just a rotated diagram of y-decay. 
The threshold energy of the reaction is roughly the same and so is also given by 



m 2 M"- 2 



1 In 



E " ,Ki \(F(r,W, ' (12) 

Just above threshold this process is also an extremely efficient method of energy 
loss. Note that while y decay destroys the incoming photon, the vacuum Cerenkov 
effect merely is an energy loss process. 

LV-modified threshold reactions: Photon Absorption 

A process related to photon decay is photon absorption, 77— > e + e~ . Unlike pho- 
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ton decay, this is allowed in Lorentz invariant QED and it plays a crucial role in 
making our universe opaque to gamma rays above tens of TeVs. 

If one of the photons has energy (Oq, the threshold for the reaction occurs in 
a head-on collision with the second photon having the momentum (equivalently 
energy) kn = m 2 / (Oq. For example, if ku = 10 TeV (the typical energy of inverse 
Compton generated photons in some active galactic nuclei) the soft photon threshold 
©o is approximately 25 meV, corresponding to a wavelength of 50 microns. 

In the presence of Lorentz violating dispersion relations the threshold for this 
process is in general altered, and the process can even be forbidden. Moreover, as 
firstly noticed by Kluzniak [85], in some cases there is an upper threshold beyond 
which the process does not occur. Physically, this means that at sufficiently high 
momentum the photon does not carry enough energy to create a pair and simul- 
taneously conserve energy and momentum. Note also, that an upper threshold can 
only be found in regions of the parameter space in which the y-decay is forbidden, 
because if a single photon is able to create a pair, then a fortiori two interacting 
photons will do [32]. 

Let us exploit the above mentioned relation rj J = (— )"T}i between the electron- 
positron coefficients, and assume that on average the initial state is unpolarized. In 
this case, using the energy-momentum conservation, the kinematics equation gov- 
erning pair production is the following [34] 

_^_.^ +t _,(^. +( . ) . (1 _, r .) 

where f = (*M/M"~ 2 and f) = 77'"' /M"~ 2 are respectively the photon's and elec- 
tron's LV coefficients divided by powers of M, < y < 1 is the fraction of momen- 
tum carried by either the electron or the positron with respect to the momentum k 
of the incoming high-energy photon and CO/, is the energy of the target photon. The 
analysis is more complicated than simple one-particle initial state decay or radia- 
tive processes. In particular it becomes necessary to sort out whether the thresholds 
are lower or upper ones and whether they occur with the equal or different pair 
momenta. 



4.2.2 Constraints from GZK secondaries 

One of the most interesting features related to the physics of Ultra-High-Energy 
Cosmic Rays (UHECRs) is the Greisen-Zatsepin-Kuzmin (GZK) cut off [45, 46], 
a suppression of the high-energy tail of the UHECR spectrum arising from inter- 
actions with CMB photons, according to py ' — > A + — > p7t°(n7t + ). This process has 
a (LI) threshold energy ~ 5 x 10 19 (afc/1.3 meV) -1 eV ((Ob is the target pho- 
ton energy). Experimentally, the presence of a suppression of the UHECR flux was 
claimed only recently [21, 20]. Although the cut off could be also due to the finite 
acceleration power of the UHECR sources, the fact that it occurs at the expected 
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energy favors the GZK explanation. The results presented in [86] seemed to further 
strengthen this hypothesis (but see further discussion below). 

Rather surprisingly, significant limits on £, = £ ( 6 ) and 77 = r\ ^ for the proton 
can be derived by considering UHE photons generated as secondary products of the 
GZK reaction[87, 88]. This can be used to further improve the constraints on dimen- 
sion 5 LV operators and provide a first robust constraint of QED with dimension 6 
CPT even LV operators. 

These UHE photons originate because the GZK process leads to the production 
of neutral pions that subsequently decay into photon pairs. These photons are mainly 
absorbed by pair production onto the CMB and radio background. Thus, the fraction 
of UHE photons in UHECRs is theoretically predicted to be less than 1% at 10 19 eV 
[89]. Several experiments imposed limits on the presence of photons in the UHECR 
spectrum. In particular, the photon fraction is less than 2.0%, 5.1%, 31% and 36% 
(95% C.L) at E = 10, 20, 40, 100 EeV respectively [90, 91]. 

However, we have just seen that pair production can be strongly affected by LV. 
In particular, the (lower) threshold energy can be slightly shifted and in general an 
upper threshold can be introduced [32]. If the upper threshold energy is lower than 
10 20 eV, then UHE photons are no longer attenuated by the CMB and can reach 
the Earth, constituting a significant fraction of the total UHECR flux and thereby 
violating experimental limits [87, 88, 92]. 

Moreover, it has been shown [88] that the y-decay process can also imply a sig- 
nificant constraint. Indeed, if some UHE photon (Ey ~ 10 19 eV) is detected by ex- 
periments (and the Pierre Auger Observatory, PAO, will be able to do so in few years 
[90]), then y-decay must be forbidden above 10 19 eV. 

In conclusion we show in Fig. 1 the overall picture of the constraints of QED 
dimension 6 LV operators, where the green dotted lines do not correspond to real 
constraints, but to the ones that will be achieved when AUGER will observe, as 
expected, some UHE photon. 



5 Dynamical Frameworks II: Rotationally Invariant CPT odd 
Standard Model Extension with mass dimension five operators 

5.1 The model 

Myers & Pospelov [38] found that there are essentially only three operators of di- 
mension five, quadratic in the fields, that can be added to the QED Lagrangian that 
give rise to dispersion modifications of type n = 3. 5 These extra-terms are the fol- 
lowing: 



5 Actually these criteria allow the addition of other (CPT even) terms, but these would not lead 
to modified dispersion relations (they can be thought of as extra, Planck suppressed, interaction 
terms) [39]. 
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Fig. 1 LV induced by dimension 6 operators. The LV parameter space is shown. The allowed re- 
gions are shaded grey. Green dotted lines represent values of (77, §) for which the /-decay thresh- 
old ky_j ec ~ 10 19 eV. Solid, blue lines indicate pairs (JJ,£) for which the pair production upper 
threshold k u „ ~ 10 20 eV. Picture taken from [93]. 



-^u'"F ma (u-D)(u n F"«) + —u'"y Y ,n(Zi + Czfc)(«--D)V, (14) 

where F is the dual of F and £12 are dimensionless parameters. All these terms 
also violate CPT symmetry. More recently, this construction has been extended to 
the whole SM [39,41,42]. 

From (14) the dispersion relations of the fields are modified as follows. For the 
photon one has 

£(3) 

col=^±"— k \ (15) 

where £ = t, and the + and — signs denote right and left circular polarisation, 
while for the fermion (with the + and — signs now denoting positive and negative 
helicity states) 

3 

Ei=p 2 + m 2 + ri2 ] {-;, (16) 

(3) 

with 77^ = 2(£i ± £2). For the anti-fermion, it can be shown by simple "hole inter- 
pretation" arguments that the same dispersion relation holds, with TJ± ' "^ = — f]^^ 
where af and / superscripts denote respectively anti-fermion and fermion coeffi- 
cients [34, 94]. Note that if CPT ends up being a fundamental symmetry of nature 
without Lorentz symmetry it would forbid all of the above mentioned CPT odd op- 
erators. 
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With the specific dimension five operators in hand, we can now return to a previously 
mentioned point, that RG flow does not significantly suppress the sizes of operators 
at low energies. Let us consider the evolution of the dimension five LV parameters. 
Bolokhov & Pospelov [39] addressed the problem of calculating the renormalization 
group equations for QED and the Standard Model extended with dimension-five 
operators that violate Lorentz Symmetry. In the framework defined above, assuming 
that no extra physics enters between the low energies at which we have modified 
dispersion relations and the Planck scale at which the full theory is defined, the 
evolution equations for the LV terms in eq. (14) that produce modifications in the 
dispersion relations, can be inferred as 

dC,\ 25 a d& 25 a 5 a e dt, la 2 a 
dt 12 % dt 12 % 12 n dt 12 % 3 % 

where a = e 2 /4n ~ 1/137 ih = 1) is the fine structure constant and t = ln(/x 2 //Xg ) 
with /i and jiQ two given energy scales. (Note that the above formulae are given 
to lowest order in powers of the electric charge, which allows one to neglect the 
running of the fine structure constant.) 

These equations show that the running is only logarithmic and therefore low 
energy constraints are robust: 0(1) parameters at the Planck scale are still 0(1) at 
lower energy. Moreover, they also show that 77 j_ and Tji cannot, in general, be 
equal at all scales. Similar calculations in the context of the renormalizable SME 
give equivalent results. 



5.3 Constraints 

We now detail some of the constraints that can be put on the CPT odd dimension 
five operators. For a thorough review of these constraints, see also [93]. 

5.3.1 Photon time of flight 

Although photon time-of-flight constraints from high energy photons propagating 
from cosmologically distant objects currently provide limits several orders of mag- 
nitude weaker than other constraints, they have been widely adopted in the astro- 
physical community. They were one of the first to be proposed in the seminal paper 
[26] . More importantly, given their purely kinematical nature, they may be applied 
to a broad class of frameworks, even beyond EFT with LV. 

In general, a photon dispersion relation in the form of (15) implies that photons 
of different colors (wave vectors k\ and £2) travel at slightly different speeds. Let us 
first ignore any birefringence effects, and just consider some coefficient <f;M that is 
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universal for all photons. Then, upon propagation on a cosmological distance d, the 
effect of energy dependence of the photon group velocity produces a time delay 

1 h n ~ 2 

At (n) = n_}_ k 2 -*i ^n) d (18) 
2 M£f 2 

which clearly increases with d and with the energy difference as long as n > 2. 
The largest systematic error affecting this method is the uncertainty about whether 
photons of different energy are produced simultaneously in the source. 

One way to alleviate systematic uncertainties is to apply EFT, which gives more 
information than just a modified dispersion for photons. In particular, one knows 
that the term generated by ^ \ in an EFT context implies birefringence. Further- 
more, photon beams generally are not circularly polarized; thus, they are a superpo- 
sition of super and subluminal circularly polarized modes. Hence one can remove 
any systematic uncertainty relating to source dynamics or measured energy by mea- 
suring the velocity difference between the two polarization states at a single energy, 
corresponding to 

At = 2\^ 3) \kd/M Pl . (19) 

This bound would require that both polarizations be observed and that no spurious 
helicity-dependent mechanism (such as, for example, propagation through a bire- 
fringent medium) affects the relative propagation of the two polarization states. 

Note that one doesn't have to actually measure the actual polarization - a double 
peak in the high energy spectra of GRB's with a separation that scaled linearly with 
distance would be a smoking gun for birefringent theories. However, if the polariza- 
tion states do not fully separate then extracting a signal becomes more complicated. 
Since current time of flight constraints compare low to high energy photons, bire- 
fringence destroys this. A GRB burst would not have low energy photons arriving 
first with high energy following (or vice vera). Instead the structure of the burst for 
a birefringent theory would be high energy photons first, low energy photons fol- 
lowing, and then more high energy photons at the end. Therefore, the net effect of 
this superposition may be to partially or completely erase the time-delay effect as it 
is usually calculated. 

In order to compute this modulating effect on a generic photon beam in a bire- 
fringent theory, let us describe a beam of light by means of the associated electric 
field, and let us assume that this beam has been generated with a Gaussian width 

E = A fe ! ^of-* + («o)z) e -(z-vJ0 2 5^g + + g i(A)'-*-(^Q)z) e -(i:-Vgf) 2 5^^_ N \ ; ( 2 fJ) 



where Qo is the wave frequency, 8Qo is the gaussian width of the wave, A: ± (i2o) 
is the "momentum" corresponding to the given frequency according to (15) and 
e± = {e\ ± iej) / y/2 are the helicity eigenstates. Note that by complex conjugation 
e\ = <?_ . Also, note that (©) = <» =p § C0 2 /M P \. Thus, 

E=Ae iti (t-z) ^ni/M M z e -(z-yp) 2 snl )i+ + e -^ni/M n z e Hz-vp) 2 sni i ^ (21) 
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The intensity of the wave beam can be computed as 

E-E* = \A\ 2 f e l%tf/Mi>iZ + e -2it;G%/Mnz\ g -5fl2(( z -v+t) 2 +(z-v-f) 2 ) 

= 2|A|V 2M o(z-0 2 cos ( 2^Q z) e- 2 ?i (5n ° t)2 . (22) 

V ™pi / 

This shows that there is an effect even on a linearly-polarised beam. The effect is a 
modulation of the wave intensity that depends quadratically on the energy and lin- 
early on the distance of propagation. In addition, for a gaussian wave packet, there is 
a shift of the packet centre, that is controlled by the square of % ' 3 ' /Mp\ and hence is 
strongly suppressed with respect to the cosinusoidal modulation. Hence by looking 
for modulation of the signal with energy and distance one can in principle determine 
if the LV is birefringent from time of flight information even if the photons arrive as 
part of one "burst". 

So far, the most robust constraints on t, ' 3 ', derived from time of flight differences, 
have been obtained within the D— brane model (discussed in section 9.1) from a 
statistical analysis applied to the arrival times of sharp features in the intensity at 
different energies from a large sample of GRBs with known redshifts [95], leading 
to limits % ( 3 ) < <9(10 3 ). The importance of systematic uncertainties can be found in 
[96], where the strongest limit /' 3 ' < 47 is found by looking at a very strong flare 
in the TeV band of the AGN Markarian 501. Finally, an extremely strong limit (for 
this method at least) of ^> < 0.8 has been obtained from the short, high energy 
gamma ray burst GRB 090510 [97]. 



5.3.2 Photon polarization 

Since electromagnetic waves with opposite circular polarizations have slightly dif- 
ferent group velocities in rotationally invariant EFT LV when CPT is violated, the 
polarization vector of a linearly polarised plane wave with energy k rotates. Dur- 
ing the wave propagation over a distance d, the rotation angle for n = 3 dispersion 
modifications is [34] 6 

8 W = ^L±Wj g ^, (24) 
2 2Mp\ 

Observations of polarized light from a distant source can then lead to a constraint 
on l^ 3 '! that, depending on the amount of available information — both on the 



6 Note that for an object located at cosmological distance (let z be its redshift), the distance d 
becomes 

d (z) = TT f 1+Z ' == di , (23) 

#o Jo ^/n A +Q,„(l+z'y 

where d(z) is not exactly the distance of the object as it includes a (1 +z) 2 factor in the integrand 
to take into account the redshift acting on the photon energies. 
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observational and on the theoretical (i.e. astrophysical source modeling) side — can 
be cast in two different ways [98]: 

1 . Because detectors have a finite energy bandwidth, Eq. (24) is never probed in real 
situations. Rather, if some net amount of polarization is measured in the band 
k\ < E < k%, an order-of-magnitude constraint arises from the fact that if the 
angle of polarization rotation (24) differed by more than 7r/2 over this band, the 
detected polarization would fluctuate sufficiently for the net signal polarization 
to be suppressed [99, 94]. From (24), this constraint is 

£0)< ^ (25) 

This constraint requires that any intrinsic polarization (at source) not be com- 
pletely washed out during signal propagation. It thus relies on the mere detection 
of a polarized signal; there is no need to consider the observed polarization de- 
gree. A more refined limit can be obtained by calculating the maximum observ- 
able polarization degree, given the maximum intrinsic value [100]: 



i7(^)=i7(O)^/(cos(20))2 3 + (sin(20))2 3 , (26) 

where iT(0) is the maximum intrinsic degree of polarization, is defined in 
eq. (24) and the average is weighted over the source spectrum and instrumental 
efficiency, represented by the normalized weight function 7(k) [99]. Conserva- 
tively, one can set 17(0) = 100%, but a lower value may be justified on the basis 
of source modeling. Using (26), one can then cast a constraint by requiring II ) 
to exceed the observed value. 

Suppose that polarized light measured in a certain energy band has a position an- 
gle o b s with respect to a fixed direction. At fixed energy, the polarization vector 
rotates by the angle (24) 7 ; if the position angle is measured by averaging over a 
certain energy range, the final net rotation (A 9) is given by the superposition of 
the polarization vectors of all the photons in that range: 

tan(2( 4 e>) = ^|i, (27) 
(cos(20))y 

where is given by (24). If the position angle at emission 0; in the same energy 
band is known from a model of the emitting source, a constraint can be set by 
imposing 

tan(2(40)) < tan(20 obs - 2Q) . (28) 

Although this limit is tighter than those based on Eqs. (25) and (26), it clearly 
hinges on assumptions about the nature of the source, which may introduce sig- 
nificant uncertainties. 



7 Faraday rotation is negligible at these energies. 
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Fig. 2 Constraint for the polarization degree. Dependence of 77 on £ for the distance of the CN 
and photons in the 150-300 keV range, for a constant instrumental sensitivity 7(k). Picture taken 
from [98]. 



In the case of the Crab Nebula, a (46 ±10)% degree of linear polarization in the 
100 keV — 1 MeV band has recently been measured by the INTEGRAL mission 
[101, 102]. This measurement uses all photons within the SPI instrument energy 
band. However the convolution of the instrumental sensitivity to polarization with 
the detected number counts as a function of energy, 7(k), is maximized and ap- 
proximately constant within a narrower energy band (150 to 300 keV) and falls 
steeply outside this range [103]. For this reason we shall, conservatively, assume 
that most polarized photons are concentrated in this band. Given afcrab = 1.9 kpc, 
ko = 300 keV and k\ = 150 keV, eq. (25) leads to the order-of-magnitude estimate 
|4| <2x 10~ 9 . A more accurate limit follows from (26). In the case of the CN 
there is a robust understanding that photons in the range of interest are produced 
via the synchrotron process, for which the maximum degree of intrinsic linear po- 
larization is about 70% (see e.g. [104]). Figure 2 illustrates the dependence of 77 on 
4 (see Eq. (26)) for the distance of the CN and for 77(0) = 70%. The requirement 
77(4 ) > 16% (taking account of a 3(7 offset from the best fit value 46%) leads to the 
constraint (at 99% CL) 

|4|<6xl0- 9 . (29) 

It is interesting to notice that X-ray polarization measurements of the CN already 
available in 1978 [105], set a constraint |4 | < 5.4 x 10~ 6 , only one order of magni- 
tude less stringent than that reported in [106]. 

Constraint (29) can be tightened by exploiting the current astrophysical under- 
standing of the source. The CN is a cloud of relativistic particles and fields powered 
by a rapidly rotating, strongly magnetized neutron star. Both the Hubble Space Tele- 
scope and the Chandra X-ray satellite have imaged the system, revealing a jet and 
torus that clearly identify the neutron star rotation axis [107]. The projection of 
this axis on the sky lies at a position angle of 124.0° ± 0.1° (measured from North 
in anti-clockwise). The neutron star itself emits pulsed radiation at its rotation fre- 
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quency of 30 Hz. In the optical band these pulses are superimposed on a fainter 
steady component with a linear polarization degree of 30% and direction precisely 
aligned with that of the rotation axis [108]. The direction of polarization measured 
by INTEGRAL-SPI in the y-rays is obs = 123° ± 11° (Iff error) from the North, 
thus also closely aligned with the jet direction and remarkably consistent with the 
optical observations. 

This compelling (theoretical and observational) evidence allows us to use eq. (28). 
Conservatively assuming 0; — o b s = 33° (i.e. 3(7 from fy, 99% CL), this translates 
into the limit 

|£ (3) | <9x 10- 10 , (30) 

and |£ (3) | < 6 x 10~ 10 for a 2a deviation (95% CL). 

Polarized light from GRBs has also been detected and given their cosmological 
distribution they could be ideal sources for improving the above mentioned con- 
straints from birefringence. Attempts in this sense were done in the past [94, 109] 
(but later on the relevant observation [110] appeared controversial) but so far we 
do not have sources for which the polarization is detected and the spectral redshift 
is precisely determined. In [111] this problem was circumvented by using indirect 
methods (the same used to use GRBs as standard candles) for the estimate of the red- 
shift. This leads to a possibly less robust but striking constraints |^ 3 '|<2.4xl0 -14 . 

Remarkably this constraint was recently further improved by using the INTE- 
GRAL/IBIS observation of the GRB 0412 19A, for which a luminosity distance of 
85 Mpc (z ~ 0.02) was derived thanks to the determination of the GRB's host galaxy. 
In this case a constraint |^' 3 '| < 1.1 x 10~ 14 was derived [1 12]. 8 



5.3.3 Synchrotron radiation 

Synchrotron emission is strongly affected by LV, however for Planck scale LV and 
observed energies, it is a relevant "window" only for CPT odd dimension five LV 
QED (and dimension four LV QED, which we describe below). We shall work out 
here the details of CPT-odd dimension five QED (« = 3) for illustrative reasons, for 
the lower dimension case see [113]. 

In both LI and LV cases [34], most of the radiation from an electron of energy E 
is emitted at a critical frequency 

3 Y 3 (E) 

co c = -eB^--^- (31) 
2 E 

where y{E) = (1 — v 2 {E))~ l l 2 , and v(E) is the electron group velocity. 



8 The same paper claims also a strong constraint on the parameter Unfortunately, such a 
claim is based on the erroneous assumption that the EFT order six operators responsible for this 
term imply opposite signs for opposite helicities of the photon. We have instead seen that the CPT 
evenness of the relevant dimension six operators imply a helicity independent dispersion relation 
for the photon (see eq.(8)). 
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However, in the LV case, and assuming specifically n = 3, the electron group 
velocity is given by 

= dE = i 3 )P_ 



Therefore, v(E) can exceed 1 if 77 > or it can be strictly less than 1 if 77 < 0. This 
introduces a fundamental difference between particles with positive or negative LV 
coefficient 77. 

If 77 is negative the group velocity of the electrons is strictly less than the (low 
energy) speed of light. This implies that, at sufficiently high energy, y{E)- <E/m e , 
for all E. As a consequence, the critical frequency CD~(y,E) is always less than a 
maximal frequency o' nax . Then, if synchrotron emission up to some maximal fre- 
quency ©obs is observed, one can deduce that the LV coefficient for the correspond- 
ing leptons cannot be more negative than the value for which o™ ax = © bs> leading 
to the bound [34] 

^>-E(^lf 2 . (33 ) 
m e \m e (o ohs J 

If 77 is instead positive the leptons can be superluminal. One can show that at 
energies E c > 8 TeV/77 1 / 3 , y(E) begins to increase faster than E/m e and reaches 
infinity at a finite energy, which corresponds to the threshold for soft VC emission. 
The critical frequency is thus larger than the LI one and the spectrum shows a char- 
acteristic bump due to the enhanced (O c . 

How the synchrotron emission processes at work in the CN would appear in a 
"LV world" has been studied in [1 14, 1 15]. There the role of LV in modifying the 
characteristics of the Fermi mechanism (which is thought to be responsible for the 
formation of the spectrum of energetic electrons in the CN [116]) and the contribu- 
tions of vacuum Cerenkov and helicity decay were investigated for n = 3 LV. This 
procedure requires fixing most of the model parameters using radio to soft X-rays 
observations, which are basically unaffected by LV. 

Given the dispersion relations (15) and (16), clearly only two configurations in 
the LV parameter space are truly different: 77+ • T]_ > and 77+ • 77_ < 0, where 7] + 
is assumed to be positive for definiteness. The configuration wherein both 7]± are 
negative is the same as the (77+ • 7j_ > 0, 77+ > 0) case, whereas that whose signs are 
scrambled is equivalent to the case (77+ • ?7_ < 0, 77+ > 0). This is because positron 
coefficients are related to electron coefficients through TJ± = — 77^ [34]. Examples 
of spectra obtained for the two different cases are shown in Fig. 3. 

A x 2 analysis has been performed to quantify the agreement between models 
and data [115]. From this analysis, one can conclude that the LV parameters for 
the leptons are both constrained, at 95% CL, to be |tj±| < 10~ 5 , as shown by the 
red vertical lines in Fig. 4. Although the best fit model is not the LI one, a careful 
statistical analysis (performed with present-day data) shows that it is statistically 
indistinguishable from the LI model at 95% CL [115]. 
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Fig. 3 Comparison between observational data, the LI model and a LV one with f] + ■ r\- < (left) 
and 7] + ■ ?7_ > (right). The values of the LV coefficients, reported in the insets, show the salient 
features of the LV modified spectra. The leptons are injected according to the best fit values p = 2.4, 
E c = 2.5 PeV. The individual contribution of each lepton population is shown. Picture taken from 
[115]. 



5.3.4 Constraints from GZK secondaries 

The same reasoning that established constraints on the CPT even higher dimension 
sector from GZK secondaries can also be applied to further strengthen the available 
constraints on CPT odd dimension 5 LV QED. In this case the absence of relevant 
UHE photon flux strengthens by (at most) two orders of magnitude the constraint on 
the photon coefficient £, ( 3 ) while an eventual detection of the expected flux of UHE 
photons would constrain rj' 3 ' for the electron/positron at the level of \r\^ \ < 10~ 16 
(see [88, 93] for further details) by limiting the gamma decay process. Note however, 
that in this case, unlike the CPT even case, one cannot exclude that only one photon 
helicity survives and hence a detailed flux reconstruction would be needed. 

5.3.5 Summary of constraints 

Constraints on LV QED at n = 3 are summarized in Fig. 4 where also the constraints 

— coming from the observations of up to 80 TeV gamma rays from the crab nebula 
[117] (which imply no gamma decay for these photons neither vacuum Cerenkov at 
least up to 80 TeV for the electrons producing them via inverse Compton scattering) 

— are plotted for completeness. 



5.4 Other threshold processes 

Once one realizes the power of the vacuum Cerenkov effect it's natural to explore 
other threshold processes that might yield useful constraints. The constraints in Fig- 
ure 4 are the best available constraints. We detail below some other processes that 
have been used to set constraints in the past. 
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Fig. 4 Summary of the con- 
straints on LV QED at order 
n = 3. The red lines are related 
to the constraints derived from 
the detection of polarized syn- 
chrotron radiation from the 
Crab nebula as discussed in 
the text. For further reference 
are also shown the constraints 
that can be derived from the 
detection of 80 TeV photons 
from the Crab nebula: the 
solid black lines symmetric 
w.r.t. the £ axis are derived 
from the absence of gamma 
decay, the dashed vertical line 
cutting the f] axis at about 
10~ 3 refers to the limit on 
the vacuum Cerenkov effect 
coming from the inferred 80 
TeV inverse Compton elec- 
trons. The dashed vertical line 
on the negative side of the 
t] axis is showing the first 
synchrotron based constraint 
derived in [114]. 



5.4.1 Helicity decay 

A slightly different version of the vacuum Cerenkov process is that of helicity de- 
cay (HD, e T — > <?±y). If ?]+ 7^ T]_, an electron can flip its helicity by emitting a 
suitably polarized photon. This reaction does not have a real threshold, but rather 
an effective one [34] — Phd = {m 2 e M / Ar]) 1 ! 3 , where Ar\ = \r]^ — 7]i 3 '| — at 
which the decay lifetime Thd changes how it scales with Arj. For Arj w 0(1) this 
effective threshold is around 10 TeV. Below threshold the lifetime, Thd is given 
by Thd > Ar]~ 3 (p/10 TeV) -8 10~ 9 s, while above threshold Thd becomes indepen- 
dent of Ar] [34] and is given by Thd ~ 10~ 10 x (10 Tev)/£ s. It is more difficult, 
however, to use helicity decay to set constraints as we do not have polarization mea- 
surements of high energy cosmic rays and the flux above a certain energy scale is 
only halved in the case of helicity decay, rather than almost completely removed as 
in the vacuum Cerenkov effect. One therefore requires more detailed knowledge of 
the source spectrum to properly apply helicity decay constraints. Nonetheless, HD 
can play a crucial role in lepton/anti-lepton propagation by basically leaving only a 
survival helicity state for each particle type. This mechanism, for example, played a 
pivotal role in the reconstruction of the Crab nebula synchrotron spectrum [115]. 




Lorentz breaking effective field theory models 



27 



5.4.2 Photon splitting and lepton pair production 

Vacuum Cerenkov radiation is an effect that would ordinarily be forbidden, however 
one can also look for modifications of normally allowed threshold reactions that are 
especially relevant for high energy astrophysics. It is rather obvious that once pho- 
ton decay and vacuum Cerenkov are allowed also the related reactions in which 
an outgoing lepton/anti-lepton pair is replaced by two or more photons, y — s- 2y 
and y — >• 3y, etc. , or the outgoing photons is replaced by an electron-positron pair, 
e~ — > e~e~e + , are also allowed. 

LV-allowed reactions: Photon splitting 

This is forbidden for |W < while it is always allowed if £jM >0 [32]. When 
allowed, the relevance of this process is simply related to its rate. The most relevant 
cases are y — > yy and y — > 3y, because processes with more photons in the final state 
are suppressed by more powers of the fine structure constant. 

The y — > yy process is forbidden in QED because of kinematics and CP conser- 
vation. In LV EFT both the kinematics changes and CP is not necessarily a good 
symmetry. However, we can argue that y^ yy is suppressed by an additional power 
of the Planck mass with respect to y — > 3y. In fact, in LI QED the matrix element is 
zero due to the exact cancellation of fermionic and anti-fermionic loops. In LV EFT 
this cancellation is not exact and the matrix element is expected to be proportional 
to at least (%E /Mp\) p , p > 0, as it is induced by LV and must vanish in the limit 

Mpi -> oo. 

Therefore we have to deal only with y— s- 3y. This process has been studied in 
[32, 118]. In particular, in [118] it was found that, if the "effective photon mass" 
m 2 , = ^Ey/M'p^ 2 <C mj, then the splitting lifetime of a photon is approximately 

T"= 3 ~ 0.025 % - 5 f- 1 (50 TeV /E y ) s, where / is a phase space factor of order 1 . 
This rate was rather higher than the one obtained via dimensional analysis in [32] 
because, due to integration of loop factors, additional dimensionless contributions 
proportional to n? e enhance the splitting rate at low energy. 

This analysis, however, does not apply for the most interesting case of ultra high 
energy photons around 10 19 eV given that at these energies m 2 ^>m 2 ifE, ^ > 10~ 17 
and <i;( 4 ' > 10~ 8 . Hence the above mentioned loop contributions are at most loga- 
rithmic, as the momentum circulating in the fermionic loop is much larger than m e . 
Moreover, in this regime the splitting rate depends only on my, the only energy scale 
present in the problem. One then expects the analysis proposed in [32] to be correct 
and the splitting time scale to be negligible at Ey ~ 10 19 eV, which therefore makes 
it not particularly competitive with other constraints. 

Lepton pair production 

The processes e~ — > e~e~e + is similar to vacuum Cerenkov radiation or helicity 
decay with the final photon replaced by an electron-positron pair. Various combi- 
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nations of helicities for the different fermions can be considered individually. If we 
choose the particularly simple case (and the only one we shall consider here) where 
all electrons have the same helicity and the positron has the opposite helicity, then 
the threshold energy will depend on only one LV parameter. In [32] the threshold 
for this reaction was derived for electron pair production, and it was found that the 
rate is a factor of ~ 2.5 times higher than that for soft vacuum Cerenkov radiation. 
Therefore, since the rate for the reaction is high as well, constraints may be imposed 
using just the value of the threshold energy. 9 



6 Dynamical frameworks III: Rotationally Invariant Minimal 
Standard Model Extension in QED 

6.1 The model 

The Minimal Standard Model Extension or mSME, is the set of renormalizable 
operators that generate LV but maintain the existing particle content of the standard 
model and do not violate gauge invariance. A subset of the mSME, the rotationally 
invariant LV operators are 

1 _ ° 1 ° 

-a« M V / 7 M V / ^^«MV / 757 Ai V / + 2 iCM M M vV / 7 M £)V W+ ^"^vWYsY^ D v y (34) 

for fermions and 

~(k F )u K r]^u v F Kl F^ + ^k AF u K e Ka p r A a F^ (35) 

for photons. The dimension 3, CPT odd Piaf term generates an instability in the 
theory and so we will set it to zero from here on out. The a term can be absorbed 
by shifting the phase of the fermion field and so we will ignore it temporarily (we 
shall return to this term when we deal with gravity). Note however that while the a 
term can be absorbed in QED, it can be measured in matter sectors where the phase 
of the fermion is important. For example, neutrinos can constrain differences in a 
between species [42]. 

The corresponding high energy (Mpi E^> m) dispersion relations for QED can 
be expressed as (see [7] and references therein for more details) 

E l\ = m l + P 2 + fe l) P + fe 2) P 2 electrons (36) 

Ey = (1 + fy 2) )P 2 Photons (37) 

9 One could of course consider any lepton/anti-lepton pair as produced, for example the reaction 
e~ — > e~vv. While standard particle physics arguments imply that the rate will be roughly equiva- 
lent to the e~e + pair production case [119] given the same order of coefficients, and the threshold 
will be slightly lower, the constraints are on a higher dimensional parameter space and so less 
useful. 
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where fg = —2bs,fg 2 ^ = — (c — ds), and /j 2 ' = kpj2 with s = ±1 the helicity state 
of the fermion [7]. The anti-fermion dispersion relation is the same as (36) with the 
replacement p — s- —p, which will change only the /j 1 ' term. 

Note that the typical energy at which new phenomenology should start to appear 
is quite low. In fact, taking for example / e ' ~ 0(1), one finds that the corresponding 
extra-term is comparable to the particle mass m precisely at p ~ m. Even worse, 
for the linear modification to the dispersion relation, we would have, in the case 
in which /j 1 ' ~ (9(1), that p th — m 2 /M P \ — 1CT 17 eV for electrons. (Notice that 
this energy corresponds by chance to the present upper limit on the photon mass, 
m r < l(T 18 eV [120].) 



6.2 Constraints 

In contrast to the higher dimension operators, LV due to the mSME does not grow 
with energy. Therefore astrophysics with its higher energies does not necessarily 
provide a tremendous advantage over laboratory and terrestrial experiments when 
testing the mSME. We list below the best constraints currently available on the QED 
sector of the mSME. For a recent listing of all the constraints as well as the non- 
rotationally invariant case, see [73]. 

6.2.1 Spin polarized torsion constraints on b for electrons 

Both spin polarized and unpolarized torsion balances can place limits on the 
mSME. Spin polarized torsion balances place limits on the electron sector of the 
mSME [121], while unpolarized torsion balances constrain the gravitational sec- 
tor [122]. Spin polarized torsion balances constrain the electron sector, as the tor- 
sion balances are constructed to have a large number of aligned electron spins (for 
example, a simple magnet attached on a torsion fiber is a very crude spin polarized 
torsion balance). Usually of course, the magnet design is optimized to search for LV. 
For example, a vertical stack of an octagonally symmetric pattern of magnets con- 
structed to have an overall spin polarization in the octagon's plane has been used as 
a torsion balance with a net electron spin polarization of 10 23 electron spins [121]. 
The mSME coefficients give rise to an interaction potential for non-relativistic elec- 
trons which produces an orientation dependent torque on the torsion balance which 
is measured using the twist of the torsion fiber. The torsion balance in its sealed 
vacuum chamber is mounted on a rotating turntable, which allows for very sensitive 
detection of any anomalous torque as a function of the rotation frequency (and the 
earth's rotation and motion in the solar system). Current limits on b from torsion 
balances are of the order of 1(T 27 GeV [121]. 
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6.2.2 Accelerator bounds on c for electrons 



Accelerator beams of various subatomic particles are produced with (roughly) time 
independent energies. Energy loss mechanisms such as the vacuum Cerenkov effect 
q — > q + y (see section 4.2. 1) for charged fermions are incompatible with a constant 
beam energy as long as the energy loss rate is high enough. Limits can therefore be 
set on the c coefficient for accelerated charged particles, as if c is too much larger 
than zero then above the n = 2 threshold energy (12) radiative energy losses be- 
come significant over the beam time. The LEP experiment accelerates electrons and 
positrons to energies of roughly 100 GeV in the lab frame. The beam energy is care- 
fully measured and it is known that the beam does not lose any significant energy to 
vacuum Cerenkov radiation. Additionally the synchrotron emission from the beam 
has been measured. Since synchrotron emission is sensitive to the Lorentz factor of 
the accelerated particles, the existence of synchrotron radiation of a certain spec- 
trum given a source particle energy and path provides constraints on any LV present 
(see section 5.3.3). Using the characteristics of the LEP beam and the associated 
synchrotron radiation spectrum limits c for electrons to be of order \c\ < 10 -15 . 

6.2.3 Astrophysical inverse Compton bounds on d for electrons 

d causes a spin dependent change to the dispersion relation for a fermion. As such, 
one of the helicity states for a fermion is always subluminal and one can't triv- 
ially apply simple astrophysical arguments based on the vacuum Cerenkov effect, 
etc. without knowing the helicity of the measured fermion. However, there is a way 
around this by considering the dynamics of sources of high energy photons pro- 
duced by inverse Compton scattering. We receive energies from the radio up to 80 
TeV from astrophysical sources such as the Crab nebula [117]. The overall spec- 
trum of these sources is well understood and the high energy emission is dominated 
by inverse Compton scattering of accelerated electrons off of lower energy photons. 
In order to constrain d one can't simply require the absence of vacuum Cerenkov 
radiation for the electrons, as half of any population will still be present under the in- 
fluence of a d term and available to inverse Compton scatter. However, the magnetic 
field present in the same sources causes the electrons to precess, thereby destroying 
any initial polarization of the electrons. Therefore, as argued in [123] both helicities 
of electrons must be present and stable if there is to be inverse Compton radiation. 
The vacuum Cerenkov effect is therefore forbidden for both helicities, which allows 
one to put double sided constraints on d of the order of \d\ < <9(10~ 12 ) [123]. 10 



10 Note that the bounds presented here are weaker by a factor of 10~ 3 , as we have used the CMB 
frame as the rest frame rather than the sun centered frame and therefore the strength of the bounds 
are weakened by the v/c of the Sun with respect to the CMB. 
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6.2.4 Cosmic ray and HESS bounds on k F 

Bounds on the fermion sector that use processes involving the fundamental QED 
vertex implicitly assume that the photon sector is unmodified. Which sector is the 
"unmodified" sector is arbitrary, in that the limiting speed of one of the sectors can 
be defined to be the "speed of light". When constraining photon coefficients, we 
simply make the opposite assumption - we assume the fermion sector is unmodi- 
fied and constrain kF- Constraints on kF can be generated for protons and photons 
using the same effects as above, i.e. the necessary absence of vacuum Cerenkov 
and photon decay if we see high energy cosmic ray protons and TeV gamma rays. 
In [124] the authors used the necessary absence of p — > p + y and a Pierre Auger 
event in conjuction with the excess of high energy TeV gamma rays observed by the 
HESS telescope (which forbids y— > p + p) to produce a two sided bound on kF of 
-9x IQ- 16 <k f <6x lfT 20 . 



7 Dynamical Frameworks IV: Gravity 
7.1 The model 

The SME is constructed by coupling matter terms to non-zero LV tensors in vacuum. 
If we left the tensors as constants, without any sort of dynamics, then one would 
break general covariance. It may be that for other reasons one may want to change 
the underlying symmetry structure for gravity in just this way; such is the case with 
Hofava-Lifshitz gravity [36]. Alternatively, if one wants to preserve general covari- 
ance one can by promoting any LV tensors to dynamical fields. Dynamical Lorentz 
breaking is also sometimes called "spontaneous Lorentz violation", although this is 
a bit of a misnomer as there are models where there really is no Lorentz invariant 
phase. There are many ways that a dynamical field can generate a LV term in vac- 
uum. If we restrict ourselves to rotational invariance, then it is natural to generate 
LV couplings by including in the action either a scalar or a timelike vector field that 
takes a vacuum expectation value. In the case of a scalar, one can use a shift sym- 
metry (0 (x) —> (x) + 0o ) to construct actions for which the derivative of the scalar 
takes a non-zero value (c.f. [125]). In the vector case, one simply puts a potential 
for the vector field such that the vector acquires a vev. We concentrate on the vector 
case here as it is the simplest model that allows for rotationally invariant Lorentz 
violation [126]. Just as the SME is a derivative expansion in derivatives of matter 
fields, one can treat the vector field in the same way. If we still denote the vector 
field by u a then we can write the low energy action for gravity and the vector as 




(38) 



where Jz?„ is given by 
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JSP, = -Z a \ 5 {V a uV){Vpu 5 ) + V{u). (39) 

The tensor Zf^yg is defined as [71] 11 

Z a \ 5 = c lg a Pg rS + c 2 8 a r 8 P s + c 3 8 a s 8\-c 4U a u Pg rS , (40) 

where c,, / = 1 , . . . ,4 are simple coefficients of the various kinetic terms and V(u) is 
a potential term that generates a non-zero vev for u a . An additional term, R a j,u a u b 
is a combination of the above terms when integrated by parts, and hence is not 
explicitly included here. 

In general, such a theory possesses four additional degrees of freedom. One of the 
vector components, however, will necessarily have a wrong sign for its kinetic term, 
thereby generating a ghost excitation. This can be remedied, at least at low energies, 
by choosing the potential V(u) = A (u 2 + 1 ), where A is a Lagrange multiplier. This 
fixes the norm of u a and removes the ghost excitation 12 so the excitations of the 
vector all have positive norm (c.f. the discussions in [70] and [127]). With this form 
of the potential, the vector theory is known as "Einstein-aether theory" [126]. We 
will refer to the vector field as an "aether" field, vector field models that couple to 
the SM Lagrangian are also sometime nicknamed "bumblebee" models. In reality 
of course the Lagrange multiplier is likely simply an approximation to a potential 
that generates dynamics such that the ghost only has UV effects. For the purposes 
of constraints, we will treat the Lagrange multiplier as "the" potential term and so 
neglect any possible ghost excitation. 

The matter couplings between aether and the standard model field content are 
the same as they were before, only now the field has dynamics so there can be posi- 
tion dependent violations of Lorentz symmetry. Interestingly, some of the couplings 
to matter that are unobservable for a single fermion field can have relevant effects 
when the aehter varies. For example, the — aw^Tj/y^y/ term in the mSME could be 
removed by making a phase change for the fermion. However, once u a is dynam- 
ical and varies with position, only a single component of the term can actually be 
removed by a phase change [128]. This leads to a new type of term which requires 
gravitational/position dependent tests in the matter sector. [128]. 

Additionally, just as constraints can be put on the coupling of matter fields to 
the LV vector field, constraints can also be put on the c, coefficients themselves by 
examining the relevant PPN parameters and comparing to known PPN limits from 
solar system and other observations. 



11 Note the indicial symmetry = Z a ^^ g . 

12 Note that one could have fixed the norm to be larger or smaller than one, however one can simply 
scale u a to have norm — 1 and absorb the scaling into the coefficients. 
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7.2 Constraints 

7.2.1 Constraints on the aether kinetic terms 

Constraints from PPN analysis 

There are two primary methods to constrain the aether kinetic terms. First, one 
can compute the "parameterized post-Newtonian" or PPN coefficients for the theory 
and compare to observational tests. The only two non-vanishing PPN parameters are 
OCi , OC2 which describe preferred frame effects. a\ and (X2 have been calculated [129] 

as 

-8(c? + cic 4 ) 

Ol = -2,2 < 41 > 

2ci — c\ + C3 

_ Oi _ (ci + 2c 3 - c 4 ) (2ci + 3c 2 + c 3 + c 4 ) 
2 (ci +c 2 + c 3 )(2-ci -c 4 ) 

Current constraints are aj < 10~ 4 and Oh < 4 x 10~ 7 [130] and so from a PPN anal- 
ysis alone there is still a large 2-d region of parameter space that remains consistent 
with available tests of GR. 

Constraints from gravity-aether wave modes 

The aether introduces three new excitations into the gravity sector (there are naively 
four, but the unit constraint removes one). These excitations strongly couple to the 
metric via the constraint. The combined aether-metric modes consist of the two 
usual transverse traceless graviton modes, a vector mode, and a scalar mode [131]. 
Gravitational wave detectors can in principle test for polarizations [132], however 
there is an additional possibility to constrain the kinetic terms. The speeds of each of 
the modes can differ from the speed of light. Hence if the speeds are less than unity, 
high energy cosmic rays will emit vacuum gravitational Cerenkov radiation [133]. 
If we denote the speeds of the spin-2, spin-1 and spin-0 modes by S2,s\ t So then we 
have [131] 

S \ = {\-c x -cif x (43) 

s \ = 2ci ~ c i+ c 3 (44) 

2(ci +c 4 )(l -c\ -c 3 ) 

s 2 = (ci+c 2 + c 3 )(2-ci-c 4 ) (45) 

(ci +c 4 )(l — c\ -c 3 )(2 + 3c2 + ci +c 3 ) 

Requiring all these speed to be greater than unity therefore puts constraints on a 
combination of the c,- coefficients. 



Combined constraints 



34 



Stefano Liberati and David Mattingly 



Even after imposing the above constraints there is still a large region of parame- 
ter space allowed. To get an estimate of the size of the space, one can set a\ and ai 
equal to zero and solve the resulting equations for C2 and C4 in terms of c\ and C3. If 
we define c + = c\ +C3 and c_ = c\ — C3 then the PPN and gravitational Cerenkov 
constraints are all satsified provided [71] 

< c+ < 1 (46) 

< c_ < C+ (47) 
3(1 -c+) 

which shows that the gravitational sector is only minimally constrained compared 
to aether-matter couplings. 



7.2.2 Gravitational constraints on aether-matter couplings 

There are couplings between aether and matter that can only be strongly tested 
when the aether is dynamical. The prime example is the coupling —aii^yY 1 y/ in the 
mSME. As said, when gravity is neglected and we're dealing with an observation 
that involves only one fermion field, then the au^ term can be absorbed by making 
a phase transformation of the fermion y/ — > e mu x yi. This cannot be done when u 
is a function of position, which would necessarily be the case if u is a dynamical 
field strongly coupled to gravity (as we've done with the constraint above). One 
can therefore use gravitational experiments to measure the position dependence of 
u 11 , similar to how redshift experiments measure the gravitational field. For a larger 
overview of constraints, see [73, 134]. We will concentrate here on a recently ap- 
plied method, that of atom interferometry [ 1 34] . 

Atom interferometry involves splitting of a beam of cold atoms such as Cs along 
two beam arms and comparing the phase difference at a point of the Cs atoms that 
travel along each beam. Since gravity and the matter sector of the mSME affect the 
phase of the atomic wavefunctions, one can look for an anomalous sidereal pattern. 
By orienting the interferometer vertically and horizontally, one can adjust the sen- 
sitivity to the gravitational and matter sectors of the mSME. Since one can also use 
other limits to control for any possible matter mSME effect on the signal one can 
thereby isolate solely gravitational effects. Hence atom interferometers have been 
used with great success on constraining the gravitational sector of the mSME. 

Consider an interferometer that is constructed to measure phase shifts of of atoms 
moving vertically vs. horizontally. There will be a contribution to the total phase 
shift A<p that comes from the gravitational redshift, i.e. some A<p rei /. In the New- 
tonian limit where is proportional to the timelike Killing vector, is given by 
uf* = (1 —U, 0,0,0) [135], where U is the Newtonian gravitational potential. If 
we couple to matter, one can easily see that the effect of the coupling via an 
au} 1 or cu^u 1 term is to change the redshift by an additional amount proportional 
to the gravitational potential, i.e. the change in frequency of a wave is given by 
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Acq = (1 + P)AU where j3 is an experimental parameter related to the fundamen- 
tal Lagrangian parameters (c.f. [136]) for a discussion. The phase shift generated 
by this extra gravitational redshift is then just A(j) rec i = (1 + j5)A(j)o, where A$o is 
the expected gravitational redshift. Constraints on /3 therefore limit the fundamen- 
tal Lagrangian parameters. /3 is limited to be zero within a few parts per billion, 
which then generates fairly strong limits on position dependence of the a and c 
couplings [134]. 



8 Neutrinos 

While we have concentrated primarily on QED, we would be remiss to not men- 
tion at least in passing the important role of neutrinos in tests of Lorentz violation. 
More so than other particles, neutrinos are uniquely suited to test various aspects 
of Lorentz symmetry. They are copiously produced in both terrestrial experiments 
and astrophysical objects, so there are solar neutrinos at < 1 MeV, controlled beams 
of neutrinos at roughly 10 GeV, atmospheric neutrinos up to 10 TeV, and (theoret- 
ically) ultra-high energy neutrinos up to 10 18 eV and above. Neutrinos are weakly 
interacting, so they propagate over long distances which allows for detailed time 
of flight and threshold reaction analyses. Neutrinos have a very small mass, which 
make threshold analyses even more sensitive. Finally, neutrinos oscillate between 
flavor eigenstates, which constrains interspecies Lorentz violation. 

The problem, of course, is that neutrinos are rather difficult to detect. However, 
thanks to the numerous ongoing neutrino experiments, many significant experimen- 
tal results have been published that can be adapted to constraining Lorentz violation 
in the neutrino sector. We now turn to the theoretical framework and current and 
future constraints on the neutrino sector. 



8.1 Neutrinos I: Species dependent Lorentz violation 
8.1.1 The model 

Neutrinos come in three distinct flavors/masses, and there is a priori no reason to 
believe that any Lorentz violating coefficients are the same for each species. It may 
even be that any Lorentz violation is not diagonal in either basis. We make a couple 
simplifying assumptions on this point - that the neutrinos are Dirac neutrinos and 
that any Lorentz violation is diagonal in the mass basis. This is a reasonable starting 
point, as one natural idea is that since any theory of quantum gravity must reduce to 
general relativity in the infrared, any Lorentz violation induced by quantum gravity 
would be primarily controlled by the charges that couple to gravity. However, this 
does not meant hat the coefficients for each mass eigenstate is the same. Indeed, one 
would expect that, due to RG effects, even if the coefficients were the same for all 
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eigenstates at one energy they would not be the same for the large range of energies 
we can test in neutrino experiments. With these assumptions, the Lorentz violating 
terms (written in the mass basis) are exactly those for the QED fermions, 

1 ° 
-a i u ll y i 'fyi - b t u^ v / ,T57 M V; + ^'^"vV;^ <5 V V* ( 48 ) 

1 ° 1 

^-idiu^u^^'f d v Wi + — -M m v7 ; 7 m (C/,i + £,275)0 ' <?) 2 y/i 

Z ZMp] 

m pi 

where ; is the mass index. We've dropped the gauge covariant derivative above, as 
it is irrelevant and couples in the flavor basis so merely would add needless com- 
plication. Also note that we have included both right and left projection operators. 
However, since standard model interactions only produce left-handed neutrinos, the 
constraints will primarily be on the corresponding left-handed operators and so we 
drop all Pr terms in our discussion of oscillations (although we will need to return 
to them when we discuss neutrino splitting). Finally, in contrast to the QED case, 
we cannot drop the a,- term as this gives a contribution to the oscillation pattern. 

The above terms and the usual Dirac Lagrangian for the neutrino yield a high 
energy neutrino dispersion relation of 

E 2 = p 2 +N 2 (49) 

Nf = m] +2( ai + b-)p - (c« + d,) P 2 + ^-af^p 2 

Mpi ' 

3 « (6 V 

The tightest constraints on species dependence comes from neutrino flavor os- 
cillation measurements. Neutrino oscillations depend on the differences in E — p 
between different neutrino eigenstates. In standard neutrino oscillations, this differ- 
ence is governed by the squared mass differences between the energy eigenstates. 
With Lorentz violation and our assumption that the Lorentz violating eigenstates 
are the mass eigenstates oscillations are governed by the differences in the effective 
mass squared. Therefore, neutrino oscillations do not probe any absolute Lorentz 
violation in the neutrino sector, but rather the differences in the dispersion relations 
between different neutrino states. 

Now let us consider a neutrino produced via a particle reaction in a definite flavor 
eigenstate / with momentum p. The amplitude for this neutrino to be in a particular 
mass eigenstate i is represented by the matrix {//,-, where £ U^Uu = 8u- The ampli- 
tude for the neutrino to be observed in another flavor eigenstate J at some distance 
L from the source, after some time T is then 

Au =Y,uy m - pL) un « L u y ,LNf /(2E)u >! ■ (5°) 
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The transition probability can then be written as 

'8NfjL\ .fSNfjL 

v^) +2C "" Si "( 



Pu = S u - £ 4F„, 7 sin 2 ( + 2G„ysin 2 -Ji- , (51) 



with SNfj = Nf-Nj and F mj and G/yy are functions of the mixing matrixes. With 
this formalism one can compare with existing neutrino oscillation experiments that 
measure Pij. Many of these experiments quote results on a deviation of the neutrino 
speed from that of light, i.e. 

-^J =E- 2 (5Nfj-Smjj) (52) 

which can be easily translated into a constraint on the coefficients in the Lagrangian, 
as we do below. 



8.1.2 Constraints 

We now turn to the constraints that can be put on the neutrino sector. For the renor- 
malizable terms governed by ai,bi,Ci,di and the non-renormalizable oc-y term that 
contributes a similar term in the dispersion relation the best known limits come from 
the miniBoone experiment [137]. They limit various combinations of coefficients at 
the 10~ 20 GeV level, which implies that the order of magnitude constraints (assum- 
ing no cancellations) on these coefficients are 

Wi\,\bi\ < l(T 20 GeV (53) 

N,K-|<io- 20 

ag < 10 9 

where we have used the mass of the neutrino as approximately 0. 1 e V and the energy 
of the neutrinos in the miniBoone experiment as approximately 1 GeV (it's actually 
slightly less). As one can see, both of the renormalizable operators are very tightly 
constrained, while the non-renormalizable operator is essentially free. 

For the higher dimension operators that contribute higher order corrections to 
the dispersion relation the best constraints to date come from the survival of atmo- 
spheric muon neutrinos observed by the former IceCube detector AMANDA-II in 
the energy range 100 GeV to 10 TeV [138]. AMANDA-II searched for a generic 
LIV in the neutrino sector [139] and achieved (4c/c)y < 2.8 x 10~ 27 at 90% con- 
fidence level assuming maximal mixing for some of the combinations Using 
the low end of the energy band (100 GeV) to be conservative, this yields order of 
magnitude constraints on the Lorentz violating coefficients of 
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ICuUCd<10- 10 (54) 

l«^l <io 7 . 

Given that IceCube does not distinguish neutrinos from antineutrinos, the same con- 
straints apply to the corresponding antiparticles. Of interest is that due to the strong 
energy dependence of the dimension six term, as more data is taken at the 10 TeV 
range the constraint will drop below 0(1) (at 10 TeV for the energy the constraint 
is already afy < 0.1). The IceCube detector is expected to improve this constraint 
to (Ac/c)ij < 9 x 10~ 28 in the next few years [140]. We also note that the lack 
of sidereal variations in the atmospheric neutrino flux also yields comparable con- 
straints on some combinations of SME parameters [141], which can be translated 
into the framework above. Finally, a nice summary of neutrino oscillation observa- 
tions, with particular attention to LIV, can be found in [142]. For a comprehensive 
listing of constraints on terms in the neutrino sector of the SME see [73, 42]. 



8.2 Neutrinos II: Species independent Lorentz violation 

Neutrino oscillation depend on 5iVy and so if all the Lorentz violating terms are 
species independent there is no contribution to the oscillation pattern from Lorentz 
violation. Hence other methods must be used to constrain these terms. The model is 
almost the same as above. The only differences are that all neutrinos have the same 
Lorentz violating coefficients and therefore the a, term can be dropped as in QED. 
Hence we proceed directly to the relevant constraints. 

8.2.1 Constraints 

Time of flight 

For pure time of flight constraints we have to date only two observations to rely 
on, the supernova SN 1987a neutrino burst and the ICARUS experiment. We deal 
first with SN1987, which was a unique event that generated the almost simul- 
taneous (within a few hours) arrival of electronic antineutrinos and photons. Al- 
though only few electronic antineutrinos at MeV energies was detected by the ex- 
periments Kamiokall, 1MB and Baksan, it was enough to establish a constraint 
(Ac/c) T0F < 10~ 8 [143] or (Ac/c) T0F < 2 x 10~ 9 [144] by looking at the dif- 
ference in arrival time between antineutrinos and optical photons over a baseline 
distance of 1 .5 x 10 s ly. Further analyses of the time structure of the neutrino signal 
strengthened this constraint down to ~ 10~ 10 [145, 146]. 

The scarcity of the detected neutrino did not allow the reconstruction of the full 
energy spectrum and of its time evolution in this sense one should probably consider 
constraints purely based on the difference in the arrival time with respect to photons 
more conservative and robust. Adopting Ac/c < 10~ 8 , the SN constraint implies the 
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following order of magnitude constraints 



\b\ < KT n GeV 
|e|,K-|<l(r 8 
a[ 5) < 10 2 1 

ICiMC 2 |<io 13 
l«i 6) | < io 34 . 



(55) 



Hence time of flight constraints are quite tight for renormalizable operators but leave 
the non-renormalizable operators effectively unconstrained. The ICARUS experi- 
ment measured the time of flight of neutrinos traveling from CERN to Gran Sasso 
in a repeat of the OPERA experiment. ICARUS found that the arrival time was 
consistent with zero and within approximately 10 ns of the expected light travel 
time [54]. The light travel time betwen CERN and Gran Sasso is roughly 2.4 ms, 
so the Ac/c from ICARUS is of the order of 10 , consistent with previous mea- 
surements made by the MINOS detector [147]. Therefore the SN1987A neutrinos 
remain the tightest constraint on Lorentz violation from time of flight experiments. 

Threshold reactions 

Threshold reactions also can be used to cast constraints on the neutrino sector. Sev- 
eral processes are of interest: neutrino Cerenkov emission v — s- y v, neutrino split- 
ting v — > wv, and neutrino electron/positron pair production V4Vc"f + . Let us 
consider for illustration the latter process as the others work similarly. Neglecting 
possible LV modification in the electron/positron sector (on which we have seen we 
have already strong constraints) the threshold energy for a dispersion modification 
that scales as p n is 



The rate of this reaction was firstly computed in [51] for n = 2 but can be easily 
generated to arbitrary n [52]. The generic energy loss time-scale then reads (drop- 
ping purely numerical factors) 



where g is the weak coupling and 0„ is Weinberg's angle. ICARUS found no 
electron-positron pair creation from the CERN neutrino beam as it passed through 
their detector [148], but the best constraint comes from the observation of upward- 
going atmospheric neutrinos up to 400 TeV by the IceCube experiment. Since the 
neutrinos propagated through the entire Earth to reach the IceCube detector the free 
path of these particles is at least longer than the Earth radius. This measurement has 




(56) 



withS (n) = £v(WM) 



n— 2 




(57) 
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been used to establish constraints on | Ci U I C2 1 < 30 while c, d are constrained at the 
10~ 12 level. No effective constraint can be optioned for the dimension six operators, 
however in this case neutrino splitting (which has the further advantage to be purely 
dependent on LV on the neutrinos sector) could be used on the "cosmogenic" neu- 
trino flux. This is supposedly created via the decay of charged pions produced by 
the aforementioned GZK effect. The neutrino splitting should modify the spectrum 
of the ultra high energy neutrinos by suppressing the flux at the highest energies and 
enhancing it at the lowest ones. In [149] it was shown that future experiments like 
ARIANNA [150] will achieve the required sensitivity to cast a constraint of order 

\ a< L ^ \ 10~ 4 . Note however, that the rate for neutrino splitting computed in [149] 
was recently recognized to be underestimated by a factor 0(E /M) 2 [119]. Hence 
the future constraints here mentioned should be recomputed and one should be able 
to strengthen the constraint by a few orders of magnitude. 



9 Other frameworks 

Specifying which dynamical framework is employed is crucial when discussing the 
phenomenology of Lorentz violations. The most robust and well-motivated frame- 
work is that which we have been discussing, effective field theory. However, it is 
not the only one and, in fact, there are reasonable arguments from holography that a 
quantum gravity theory should not necessarily be a local field theory in the UV (c.f. 
the discussion in [151]). Hence Lorentz violation may enter into low energy physics 
in novel ways. In addition, if one believes that the fundamental quantum gravity the- 
ory should not be Lorentz invariant, then one might look for ways in which Lorentz 
violation might appear outside the realm of effective field theory and so avoid many 
of the constraints that exist in the EFT framework. For completeness, and because 
the EFT approach is nothing more than a highly reasonable, but rather arbitrary "as- 
sumption", it is worth studying and constraining additional models, given that they 
may evade the majority of the constraints discussed in this review. 



9.1 D-brane models 

Ellis et. al. developed a model [56, 58] in which modified dispersion relations derive 
from the Liouville string approach to quantum space-time [152]. Liouville-string 
models [152] motivate corrections to the usual relativistic dispersion relations that 
are first order in the particle energies and that correspond to a vacuum refractive 
index 77 = 1 — (E/Mp\) a , where a = 1. Models with quadratic dependences of the 
vacuum refractive index on energy: a = 2 have also been considered [63]. 

In particular, the D-particle realization of the Liouville string approach predicts 
that only gauge bosons such as photons, and not charged matter particles such 
as electrons, might have QG-modified dispersion relations. This difference occurs 
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since excitations which are charged under the gauge group are represented by open 
strings with their ends attached to the D-brane [153], and that only neutral excita- 
tions are allowed to propagate in the bulk space transverse to the brane [154]. Thus, 
if we consider photons and electrons, in this model the parameter 77 is forced to 
be null, whereas | is free to vary. Even more importantly, the theory is CPT even, 
implying that vacuum is not birefringent for photons (| + = 



9.2 New Relativity Theories 

Lorentz invariance of physical laws relies on only a few assumptions: the principle 
of relativity, stating the equivalence of physical laws for non-accelerated observers, 
isotropy (no preferred direction) and homogeneity (no preferred location) of space- 
time, and a notion of precausality, requiring that the time ordering of co-local events 
in one reference frame be preserved [155, 156, 157, 158, 159, 160, 161, 162]. 

All the realizations of LV we have discussed so far explicitly violate the principle 
of relativity by introducing a preferred reference frame. This may seem a high price 
to pay to include QG effects in low energy physics. For this reason, it is worth 
exploring an alternative possibility that keeps the relativity principle but that relaxes 
one or more of the above postulates. 

For example, relaxing the space isotropy postulate leads to the so-called Very 
Special Relativity framework [163], which was later on understood to be described 
by a Finslerian-type geometry [164, 165, 166]. In this example, however, the gen- 
erators of the new relativity group number fewer than the usual ten associated with 
Poincare invariance. Specifically, there is an explicit breaking of the 0(3) group 
associated with rotational invariance. 

One may wonder whether there exist alternative relativity groups with the same 
number of generators as special relativity. Currently, we know of no such generaliza- 
tion in coordinate space. However, it has been suggested that, at least in momentum 
space, such a generalization is possible, and it was termed "doubly" or "deformed" 
(to stress the fact that it still has 10 generators) special relativity (DSR). Even though 
DSR aims at consistently including dynamics, a complete formulation capable of 
doing so is still missing, and present attempts face major problems. Thus, at present 
DSR is only a kinematic idea. 

Finally we cannot omit the recent development of what one could perhaps con- 
sider a spin-off of DSR that is Relative Locality, which is based on the idea that 
the invariant arena for classical physics is a curved momentum space rather than 
spacetime (the latter being a derived concept) [167]. 

DSR and Relative Locality are still a subject of active research and debate (see 
e.g. [168, 169, 170, 171]); nonetheless, they have not yet attained the level of matu- 
rity needed to cast robust constraints 13 . 



Note however, that some knowledge of DSR phenomenology can be obtained by considering 
that, as in Special Relativity, any phenomenon that implies the existence of a preferred reference 
frame is forbidden. Thus, the detection of such a phenomenon would imply the falsification of 
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10 Discussion and Perspectives 

As we have seen, for rotationally invariant QED Lorentz symmetry is extremely well 
tested, with strong constraints all the way up to dimension six operators. Lest the 
reader get a false impression, we also should mention that many other sectors of the 
standard model, from neutrons to mesons both with and without rotation invariance 
also have tight limits set on any possible deviation from Lorentz invariance. 

There are two areas where there remains immediate work to be done. First, there 
is one current caveat in regards to the dimension six operator constraints that needs 
to be resolved. As we have seen, the dimension six constraints mostly rely on the 
physics of the GZK feature of the UHECR spectrum. More specifically, UHECR 
constraints rest upon the hypothesis, not in contrast with any previous experimental 
evidence, that protons constitute the majority of UHECRs above 10 19 eV. Recent 
PAO [172] and Yakutsk [173] observations, however, showed hints of an increase 
of the average mass composition with rising energies up to E ps 10 19,6 eV, although 
still with large uncertainties mainly due to the proton-air cross-section at ultra high 
energies. Hence, experimental data suggests that heavy nuclei can possibly account 
for a substantial fraction of UHECR arriving on Earth. Furthermore the evidence 
for correlations between UEHCR events and their potential extragalactic sources 
[86] — such as active galactic nuclei (mainly blazars) — has not improved with 
increasing statistics. This might be interpreted as a further hint that a relevant part 
of the flux at very high energies should be accounted for by heavy ions (mainly iron) 
which are much more deviated by the extra and inter galactic magnetic fields due to 
their larger charge with respect to protons (an effect partially compensated by their 
shorter mean free path at very high energies). If consequently one conservatively 
decides to momentarily suspend his/her judgment about the evidence for a GZK 
feature, then he/she would lose the constraints at n = 4 on the QED sector 14 as well 
as very much weaken the constraints on the hadronic one. 

Assuming that current hints for a heavy composition at energies E ~ 10 196 eV 
[172] may be confirmed in the future, that some UHECR is observed up to E ~ 
10 20 eV [175], and that the energy and momentum of the nucleus are the sum of 
energies and momenta of its constituents (so that the parameter in the modified dis- 
persion relation of the nuclei is the same of the elementary nucleons) one could 
place a first constraint on the absence of spontaneous decay for nuclei which could 
not spontaneously decay without LV. 15 Such a constraint would place bounds on 



both special and doubly-special relativity. An example of such a process is the decay of a massless 
particle. 

14 This is a somewhat harsh statement given that it was shown in [174] that a substantial (albeit 
reduced) high energy gamma ray flux is still expected also in the case of mixed composition, so 
that in principle the previously discussed line of reasoning based on the absence of upper threshold 
for UHE gamma rays might still work. 

15 UHE nuclei suffer mainly from photodisintegration losses as they propagate in the intergalactic 
medium. Because photodisintegration is indeed a threshold process, it can be strongly affected by 
LV. According to [176], and in the same way as for the proton case, the mean free paths of UHE 
nuclei are modified by LV in such a way that the final UHECR spectra after propagation can show 
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subluminal protons, because in this case the energy of the emitted nucleon is low- 
ered with respect to the LI case until it "compensates" the binding energy of the 
nucleons in the initial nucleus in the energy-momentum conservation. A comple- 
mentary constraint, an upper limit on superluminal protons, can be obtained from 
the absence of vacuum Cerenkov emission. If UHECR are mainly iron at the high- 
est energies the constraint is given by C P < 2 x 10 2 for nuclei observed at 10 19 6 eV 
(and Cp < 4 for 10 20 eV), while for helium it is < 4 x 1(T 3 [176]. 

Second, the gravitational sector of Lorentz symmetry violation is currently less 
constrained. A large region of the parameter space of Einstein-aether theory re- 
mains unconstrained, while atom interferometry is starting to probe matter-tensor 
couplings that are sensitive to the dynamics of the Lorentz violating tensors. Useful 
modified gravity theories that also evade solar system tests, such as galileons [177], 
or allow for better UV behavior, such as the aforementioned Hofava-Lifshitz grav- 
ity, yield interesting Lorentz symmetry violating gravitational phenomenology. 
Hence fully exploring Lorentz violation in the gravitational sector is currently an 
important area that requires further progress. 

Acknowledgements We wish to thank Luca Maccione for useful insights, discussions and feed- 
back on the manuscript preparation. 
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